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We present here some studies on noise-induced order and synchronous firing in a system of bidirectionally 
coupled generic type-I neurons. We find that transitions from unsynchronized to completely synchronized 
states occur beyond a critical value of noise strength that has a clear functional dependence on neuronal 
coupling strength and input values. For an inhibitory-excitatory (IE) synaptic coupling, the approach to a 
partially synchronized state is shown to vary qualitatively depending on whether the input is less or more 
than a critical value. We find that introduction of noise can cause a delay in the bifurcation of the firing 
pattern of the excitatory neuron for IE coupling. 
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1. INTRODUCTION 

Among the host of interesting dynamical behaviour 
observed in coupled nonlinear systems, synchronization 
phenomena are one of the most fascinating and among 
the most studied. These phenomena abound in nature 
and in daily life: studies have been made of synchro¬ 
nization in numerous systems in pendulum clocks, 
electronic circuits, chemical systems, swarms of light- 
emitting fireflies, biological rhythms, coupled Josephson 
junctions, cardiorespiratory interactions, neuronal en¬ 
sembles in sensory systems, etc. to name a few. Syn¬ 
chronous activity in neuronal ensembles in particular, 
have received a great deal of attention and numerous 
studies have been done on the spiking patterns of neu¬ 
rons and neuronal networks^:! 7 . Neuronal ensembles per¬ 
form complex tasks and can extract different kinds of 
features from the information they receive and it is be¬ 
lieved that cognitive tasks such as feature extraction and 
recognition, sensory perception, etc., are brought about 
by synchronous neuronal activity. The mechanisms and 
dynamics of synchronization phenomena in neurons are 
therefore of great interest^- 1 '. 
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Neurons are of many kinds (see for instance^), but they 
may be broadly classified into two types based on their 
excitability mechanisms. This classification of neurons 
on the basis of the firing patterns of the axons was first 
done by Hodgkin & Huxley in their classic worbi. 

Neurons of type-II begin firing at a relatively high 
frequency through a subcritical Hopf bifurcation. 

Type-I neurons however can fire at arbitrarily low 
frequencies at low values of the applied input current, 
and the frequency increases with increasing values of the 
input. The equations which describe type-I dynamics 
can be reduced to the canonical normal form for a 
saddle-node bifurcation^. In appropriate parameter 
regimes the system operates in the close vicinity of a 
saddle-node bifurcation on an invariant circle and this 
ensures repetitive firing. It is interesting to study type-I 
neurons as almost all mammalian neurons fall under this 
class of excitability. 

A variety of synchronous phenomena occur in nonlin¬ 
ear systems^, but among these, complete synchroniza¬ 
tion (henceforth referred to as CS) is perhaps one of the 
most interesting^! as it brings about synchronization in 
the phases, frequencies and amplitudes of the signals. 
Although CS is widely believed to happen between iden¬ 
tical systems, and CS between non-identical systems is 
not usually expected, we have shown in a different worbi^ 
that in an ensemble of two hundred non-identical bidirec¬ 
tionally coupled type-I neurons with random strengths 
coupled in an all-to-all way, near complete synchroniza¬ 
tion does indeed unexpectedly occur under the influence 
of weak Gaussian white noise. 

In fact, a more extensive study of noise-induced 
synchronous activity in a system of just two bidirec¬ 
tionally coupled identical (excitatory-excitatory) type-I 
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neurons reveals that the noise-induced CS state gives 
way beyond a critical coupling and noise strength, to a 
desynchronized regime and then to subsequent locking 
to a different partially synchronized state. We find that 
the value of the critical coupling constant is proportional 
to the square root of the noise strength^. We were 
unable to achieve CS (noise-induced or otherwise) in a 
system of only two bidirectionally coupled non-identical 
inhibitory-excitatory neurons, although we showed in^ 
that in a larger ensemble, non-identical neurons can 
exhibit noise-induced CS. 

In their studies of type-I neurons coupled through 
time-dependent synapses, Borgers and Kopell 11 ’ 12 dis¬ 
cussed the effects of random connectivity on synchroniza¬ 
tion and the PING mechanism in networks of excitatory 
(E) and inhibitory (I) neurons. 

In this work, we report some interesting observations 
we have made in some computer studies of coupled 
generic type-I neurons when subject to weak additive 
Gaussian white noise. The dynamics of the synapses cou¬ 
pling the neurons are governed by ordinary differential 
equations which depend upon the outputs of the presy- 
naptic neurons- 11,12 — we have considered both EE and 
IE bidirectional couplings. 

We find that the coupled system shows type-I firing pat¬ 
tern for EE coupling which suggests the existence of a 
contraction region close to the stable manifold of the sad¬ 
dle — CS can thus be expected in EE coupled systems. 
This firing pattern is not seen for all the coupled neurons 
in IE systems — the existence of a contraction region is 
therefore in question in this case, giving a good reason 
not to expect CS for IE couplings. 

As pointed out in — , noise-assisted CS occurs in 
EE (identical) systems in certain parameter regimes 
beyond which CS is destroyed but replaced by a different 
partially synchronized state. 

For such a system, we obtain here the functional depen¬ 
dence between the critical value of the noise strength 
and the strength of synaptic coupling and the externally 
applied input at the point when desynchronization gives 
way to CS. 

For the system of two coupled IE neurons, noise-assisted 
CS could not be achieved, but we find that noise induces 
locking to a partially synchronized state at higher 
coupling strengths and external input values. For the 
IE system, there exists a transition point at a critical 
external-input value before which desynchronization 
increases with increasing noise strength, and beyond 
which desynchronization decreases with increasing noise 
strength. 

We also report here noise-induced delay in the bifur¬ 
cation point associated with the firing pattern of the 
excitatory neuron in the IE system. We suggest that the 
fast-slow relaxation dynamics of the neuronal inputs in 
IE systems make CS of the outputs hard to achieve. 

We organize our presentation in the following way. In 


Section 2 we review the mathematical description of the 
mechanism of firing in type-I neurons and the framework 
for studying them when they are coupled. In Sections 
3, 4 and 5 we describe our study of the coupled system 
in the presence of (weak) Gaussian white noise for EE 
and IE bidirectional couplings (Sections 4 and 5, respec¬ 
tively), separately. 

For each type of coupling we make different studies to 
get insights upon the mechanisms governing synchroniza¬ 
tion. These are discussed in separate subsections: for an 
intuitive understanding of the dynamics underlying the 
approach to synchrony or order, in Section 4.1 we draw 
an analogy with the mechanical system of two damped 
nonlinear spring-mass systems coupled together. In Sec¬ 
tion 4.2, we describe our study of CS using plots between 
instantaneous frequencies and instanataneous phases of 
the inputs to the neurons and explain why these can be 
used as indicators of CS in coupled systems. In Section 
4.3 we show how firing frequency-input plots for the cou¬ 
pled system can give indications of whether or not the 
coupled system is likely to show CS. The largest and 
second largest Lyapunov exponents have been calculated 
for the deterministic coupled system. In Section 4.4, we 
present our study of noise-induced variation of firing fre¬ 
quency. In the IE case, we show in the corresponding 
Subsection, the very interesting finding of noise-induced 
delay in the bifurcation associated with firing for the ex¬ 
citatory neuron. In Section 4.5, we present our study of 
variation of the synchronization error for CS, with noise 
strength, input, and coupling strength. We obtain useful 
functional relations between these for EE as well as for 
IE couplings. Results of similar studies are presented for 
IE coupling in Section 5 ; in Section 5.4 we report our 
very interesting observation of noise-induced delay of bi¬ 
furcation associated with the neuronal firing. In Section 
6, we show for EE coupling how the synchronization time 
varies with noise strength and suggest the possibility of 
exploiting the results for feature extraction. We conclude 
by summarizing the results of our work in Section 7. 


2. COUPLED THETA NEURONS 

In a type-I neuron, the mechanism of excitability comes 
about as follows^,—. For low values of the control pa¬ 
rameter (the input current), there is a stable fixed point 
and an unstable fixed point which is a saddle. Trajecto¬ 
ries forming the unstable manifold leave the saddle point 
and enter the stable fixed point, forming a loop in phase 
space. This closed loop therefore now contains the stable 
and unstable fixed points. As the value of the control 
parameter is increased, the two fixed points come closer 
together, and at the critical value, they coalesce and then 
disappear, resulting in a stable periodic solution corre¬ 
sponding to repetitive firing. 

The activity Xi of the ?’th type-I neuron in an ensemble 
of N neurons, can be related to the membrane conduc- 
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tance and its dynamics can be described by: 

N 

Xi = X T /j, Ii — (3i T 3 ' THj9jiSji 5 (1) 

j =1 

where we have chosen to work in units in which the time 
constant for the membrane potential is set to unityi^. 
Its total input Ii comprises of the external input /3j 
(which in this work we have taken to be constant), and 
the contributions from the presynaptic neurons, with 
i = 1,... N. The synaptic gating variable Sij represents 
the fraction of ion channels open in the jth presynap¬ 
tic neuron, and its values always lie in the range 0 to 
1, reaching the maximal value when the neuron spikes. 
aj = +1 denotes an excitatory synapse and aj = — 1 
models an inhibitory synapse, while gji denotes the mea¬ 
sure of the synaptic strength from neuron j to neuron j; 
in our work we have considered gu = 0. Eqn.(l) has no 
fixed points for p > 0. Since any solution of the equation 
tends to infinity in a finite time, a nonlinear transforma¬ 
tion to new variables 9j may be madei^: Xi = tan di/2 
which maps the real line onto a circle, and so avoids this 
blow-up. In terms of the new variables, eqn.(l) becomes 


ddj 

dt 


N 

(1 - cos 9i) + {/3i + '^2 a j9jiSji + cos 6 l ) (2) 

3 =1 


The point 9 = tt represents the point at infinity in eqn.(l) 
and is interpreted as firing of a spike; this transformation 
gives the theta neuron its name. 

The dynamics of the synapse .s J7 ; (t) follow from the dif¬ 
ferential equation considered iniir— : 

dSjj _ Sjj | c -r7(l+cosfl,) 1 ~ S ji ^ 

dt Tji T r 

The synaptic decay time has been denoted by Tji and tr 
denotes the synaptic rise time. 

It is clear from eqn.(2) that the activity of the ith neuron 
is regulated by feedback from ,s 7 ,; because Sji depends 
upon 9j which in turn depends upon 9i (j y i) through 
the coupling term having Sjj(9i). The coupled system 
under study is depicted in Fig.(l). 

Since the time dependence of the bifurcation parameter 
Ii is brought about because of the feedback a j9ji s ji 

from the other neurons received through the synapse, we 
chose in this study to keep the external input /3j constant 
throughout so that the dynamics governing the feedback 
could be understood better. 


3. COUPLED THETA NEURONS IN THE 
PRESENCE OF WEAK NOISE 

At the cellular scale, the activity of real neurons can 
be influenced by a number of other factors : physical 
as well as biochemical, and also by random influences: 


for instance, fluctuations in the neurotransmitter lev¬ 
els at synapses, conductance fluctuations in ion chan¬ 
nels or thermal fluctuations^. It is therefore of great 
importance to study the influence of random noise on 
coupled neuronal dynamics. In this work we have used 
Gaussian white noise £(f) to model these random influ¬ 
ences on the dynamics of the coupled generic neurons in 
eqn.(2). We consider neurons which are coupled bidi¬ 
rectionally as shown in Fig.(l) and subject to Gaussian 
white noise £(t) with the following properties: (£(f)) = 0 
, (£(t)£(^)) = 2 aS(t — t'), The stochastic variables are 
taken to obey Stratonovich calculus. Addition of Gaus¬ 
sian white noise £ to eqn.(l), manifests as multiplicative 
noise in eqn.(2) because of the change of variables to 9, 
so that the equations now take the form 


d9j 

dt 


N 

(1-cos 9j ) + (/3j+aj gSjj {9j )+£ (t )) (1+cos 9 t ) 

3 = 1 


( 4 ) 

We consider the neuronal output to be described by the 
variable u* = (1 — cos#i)/2 (as inH) since its time series 
resembles that of the membrane potential in real neu¬ 
rons. This transformation maps the resting point Xi = 0 
corresponding to 9i = 0 to tq = 0, and the spiking point 
9i = tt to m = 1 via the relation m = xf/(l + x ?). This 
choice of variables enables us to get some new insights 
into the dynamics underlying the onset of complete syn¬ 
chronization. Eqns.(4) and (3) then become 


tli 2(uj -\- (/3^ -|- 3 ^ T £)(1 ^3^ \!‘U’i (1 

j =i 

Sji = + exp (-2r;(l - Uj))— —(6) 

Tji Tr 

Among the different kinds of synchronous phenom¬ 
ena, we find CS most interesting because the trajectories 
of the mutually synchronizing systems become identical. 
The literature contains several papers on CS in different 
coupled systems; however systems with more complicated 
couplings, such as that described by eqns.(5) and (6) have 
not been intensively studied yet — an adequately satis¬ 
fying explanation of conditions and circumstances under 
which CS can occur in such systems is still lacking. 
Noise-induced CS in systems through common additive 
white noise was studied in?4 who showed that the ex¬ 
istence of a significant contraction region in phase space 
was a necessary condition for occurrence of CS — the sys¬ 
tems studied in their work were the Lorenz and Rossler 
systems — far easier to handle analytically in comparison 
to our system. 

For type-I neurons, we could expect CS among neurons 
through common noise alone, because of the existence of 
a contraction region close to the stable manifold of the 
saddle. When any two type-I neurons are coupled to¬ 
gether as in eqns.(5),(6), then the nature of the eigenval¬ 
ues of the stability matrix of the coupled system at the 
fixed points would determine whether such a contraction 
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region could exist, and if so, its nature. 

However, for our system in eqns.(5),(6), linear stability 
analysis fails as the Jacobian becomes singular at the 
fixed points. Physically, this is a reflection of neuronal 
spiking with the inbuilt phase-resetting mechanism con¬ 
tained in eqns.(3). Moreover, the time dependence of the 
control parameter and its continuous dependence upon 
the feedback also makes an analytical approach non- 
straightforward. We therefore perform some numerical 
and computer-based studies which provide some useful 
insights on the coupled dynamics. 

Just as in2i, we define CS to take place between neu¬ 
rons 1 and 2 on obtaining a vanishing value for the quan¬ 
tity (|Mi — W 2 I) which is the synchronization error aver¬ 
aged over all iterations. 

It is well known that inhibitory couplings help in 
producing phase-synchronous activity in deterministic 
systems^— . For the system under study in Fig.(l) too, 
we found that antiphase states are stable in the noise¬ 
less situation for EE synapses: these states become com¬ 
pletely in-phase in the presence of noise. This does not 
happen in general for IE synapse s 10 : 15 . 

In the uncoupled system neuronal firing starts after 
the threshold value /3 = 0 is crossed. In the coupled 
system under study (Fig.(l)), the picture changes and 
the input p acquires time dependance; neuronal firing is 
then controlled by the different parameters in eqns.(2) 
& (3): the coupling strength a^gji , the decay time of 
the synapse r^, the external input Pi, etc.. 

We proceed as follows in our investigation. We look at 
the behaviour of the neuronal output variables u \, U 2 for 
different parameters in the system in the noiseless case. 
The addition of noise changes the behaviour, introducing 
order into the system. However, the outcome is different 
for the EE and IE cases. To avoid any confusion, we 
therefore detail our results separately, for EE and IE, re¬ 
spectively, in the coming subsections. A fruitful way of 
understanding the physics of any system is to make an 
analogy with a more common, familiar one. Rewriting 
equations (5) and (6) enables a resemblance of the cou¬ 
pled neurons to a nonlinear anharmonic oscillator to be 
observed, to help understand its behaviour. We do this 
below. 


4. EE SYNAPTIC COUPLING 

In Fig. (2), we depict for EE coupling, the variations of 
both the neuronal output variables u\ and U 2 with the 
different parameters in the system for the noiseless case: 
the constant inputs Pi = /? 2 , the coupling strengths 
912 = 9211 and the synaptic decay time , varying one at 
a time and keeping the other two fixed. These portrayals 
of the system’s behaviour with changing parameters 
have been done differently, with the colour coding 
showing the varying values of the control parameter. 
The advantage of plotting them this way is essentially to 


portray the maximum possible information (for instance 
variation of both neuronal variables simultaneously) 
at a glance. The diagram distinguishes between states 
where both neurons are active and evolving (“firing” in 
this context implying that there is a difference between 
consecutive outputs for both neuron 1 and neuron 2), 
and states where one or both neurons are quiescent or 
have reached a steady (unchanging) state (“non-firing”) 
implying that there is no difference between consecutive 
outputs of one or both neurons). Thus the following 
states would be classified as “firing” in the plot: 

(i) | Ui(t + 1) — Ui(t)\ > 0 and | Uj(t + 1) — Uj(t)\ > 0. 

* 7^ j i 

(ii) | Ui{t + 1) — Ui(t )| = 0 and | Uj(t + 1) — Uj(t)\ = 0, 
*7 ^ 3 ■ 

Thus, even the state where u\(t + 1) = 1, ui(t) = 1, 
with |it 2 (t + 1) — U 2 (t)\ > 0, i.e., where neuron 1 has 
reached its peak value of 1 in consecutive firings, while 
the other neuron may have likewise peaked or have 
shown some or no change, we would classify it under the 
“non-firing” plot. Thus, it would be most appropriate to 
term these plots “Activity Diagrams”, indicative of the 
state of active evolution for both neurons in the coupled 
system. 

In the noiseless case the neurons can be excitable and 
fire only if the inputs Pi are above the threshold Pi = 0. 
For Pi < 0, there are stable and unstable fixed points, 
but no periodic firing can arise as they are below the 
threshold. The maximal value of unity for ix, corresponds 
to the generation of the spike (peak value). 

In Fig.(2a), we observe in the first plot where Pi = P 2 
is the varying parameter for fixed gij and ly,-, that above 
the firing threshold p > 0, the points in the plot for any 
P trace out a quadrant of a circle — this is especially 
apparent for large p, (P ss 4) : this is a reflection of 
the stability of the antiphase states u\ and 112 in the EE 
case. In the adjacent (middle) plot where gij = g rl is the 
control parameter, all the points plotted are the square 
“firing” ones because the input has been chosen in this 
case to be Pi = P 2 =0.1 > 0. Here again it was seen 
that the solutions ui and 112 are out of phase and the 
points trace out a roughly circular curve: The last plot 
in Fig.(2a) where the decay time Tij = Tji of the gating 
variable is the parameter shows that the magnitudes of 
Mi and M 2 do not show appreciable change with change 
in magnitude of T t j. 

To make the plots clearer, we have plotted the variation 
of both neurons zq with p in the standard manner, in 
Fig. (2b). For p < 0 the neurons do not fire and each 
solution is non-oscillatory; after the threshold is crossed 
however, i.e., for P > 0, both neurons show oscillatory 
behaviour — the plot therefore shows bifurcation of neu¬ 
ral activity. We observe that Mi and M 2 oscillate out of 
phase with each other. 

In Fig. (3) we show the development of noise-induced 
synchrony as the noise strength is increased from 0.375 
(left) to 1.0 (right) for EE coupling: a drastic change in 
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the firing activity pattern is observed on introduction of 
weak noise: in contrast to Fig.(2), all neurons now fire. 
Even for negative values of /3i & (corresponding to 
stable points in the flow of the uncoupled system: the 
neurons do not fire in the deterministic system), addi¬ 
tion of weak noise induces firing by raising total input 
levels to threshold values. One observes that the nearly 
symmetrical distribution of points on either side of the 
diagonal u\ — u 2 = 0 in (Figs.(2a)) in the noiseless case, 
converges towards the diagonal in the presence of noise. 
As the noise strength is increased, the firing of the two 
neurons is completely synchronized and for a = 1.0, all 
the points line up along the diagonal and are CS solu¬ 
tions. 


lead to both neurons eventually firing at the same rate if 
the neurons happen to be identical. 

The addition of common weak noise £(f) to both neu¬ 
rons (eqns.(5)) changes the potential at the synapses (in 
mechanistic terms, it affects the “stiffness” of the synap¬ 
tic conductances in eqns.(7)), and also modulates the fir¬ 
ing thresholds of u± and u 2 by injecting energy at random 
points in time. 

For EE coupling, the thresholds of each neuron is modu¬ 
lated similar to the other due to the common noise and if 
the noise strength is adequate, neuronal firing eventually 
occurs in synchrony. 

4.2 Variation of instantaneous frequency with 
instantaneous phase of neuronal inputs 


4.1 Analogy to a mechanical system 

We can get some insight into the dynamics underlying 
the development of CS in the EE systems by drawing an 
analogy with a mechanical system. To do this, we rewrite 
eqns.(5),(6) for N = 2, aq = 1, a 2 = 1 as: 

Ul = !d{ui) + filfniui) + g 2 l(s 2 l + l)(li2)/n(wi) 

U2 = fD(u 2 ) + p2 fn(u 2 ) + ffl 2 (Sl 2 + 1 ) )f n (u 2 ) 

S21 + h{u 2 )s21 + ($(u 2 ) + ^(W2)512(S12 + 1))S21 

= -'f'(u 2 )S2l£,(t) 

S12 + h(ui)s 12 + (<f>(ui) + ^(wi)52l(S21 + 1))S12 

= -^(Ul)si2^(f) 

where 


To understand how this CS actually comes about, 
ini£ we constructed the analytical signal Bi(t)e’ lpi ^ = 
Jj(t) + iH(Ii(t)) for the inputs It and those for the neu¬ 
ron outputs Ui : Wi(t) = Ui(t ) + H(m(t )) = Ri(t)e l<l,i W 
using Hilbert transforms and showed that the instanta¬ 
neous amplitudes Bi and instantaneous phases pi of the 
S,fn(u\) i n p U ts evolve in time according to: 

Un(u 2 ) 


Biit) — —Bi( -1- 

V Tii 


0 —2rj(l—Rj cos <fij) 


(7) 


3 1 

A 

-COS Pi 

T ji 

1 


TR 


—(a 

tr V 


+ 


\ " a j9ji 


tr 


j 


cos 


/n('«i) = 2zq /2 (l - Ui) 1/2 ; f n {ui) = 2rq /2 (l - w l ) 3/2 ; 
1 e —2rj(l— Ui) , 

- + 

'Ti, 


x e 


h(ui) = + 

$( ui ) = 2 77 


u t r 

e -2r,(l- Ui ) 


Pi{t) = — 


—277(1 — Rj cos (f>j) 
0 —2rj(l—Rj cos <f>j) 


■ cos(2r]Rj sin < p 3 


((2?? Rj sin <f>j) - p^j 


A 


tr 


■ s\\\(2riRj sin (f>j ) — sin pi 
■fcsi 7 ~j i 


= 2rj- 


tr 

=—297(1—-Ui) 


(/D(Mi) + A fn(Ui)) 


+ — g: (a + sin sin A) - Pi) 


Tr 


-fn{Ui) 


( 8 ) 


x e 


-2rj(l — Rj cos cf)j ) 


(9) 


and where a shift of by unity has been performed 
in eqns.(7) (and later on in eqns.(14)) for the sake of 
putting these in the more convenient form displayed. As 
can be seen, these equations are not easily amenable to 
a quick analysis of the physics that they describe. The 
equations for the gating variables Sij in eqns.(7) resemble 
those of a nonlinear anharmonic oscillator with a nonlin¬ 
ear damping term and a nonlinear restoring force term 
with variable spring coefficient ($ (m) ±(ui )gji (sji +1)), 
with the plus sign for excitatory couplings and minus 
sign for inhibitory couplings. In the EE system, the vari¬ 
able spring coefficient is (<$>( 7 ^) + ^(u^gjiSji) for both 
S21 and S12 and it is like a system of two damped (non¬ 
linear) spring-mass systems coupled together, with both 
“springs” getting extended or both getting compressed 
(though by different amounts) at the same time. This, 
coupled with the same (plus) sign in the Ui equations 


We argued irai£ that CS is achieved if the rate of change 
of instantaneous frequency (3) with the instantaneous 
phase (<f>) of the input for neuron 1 matches exactly that 
for neuron 2. This is not a trivial statement since the 
system incorporates feedback as well as noise. 

This matching of values happens because noise produces 
delay in the decay time of the synaptic gating variables, 
causing CS to eventually occur. In Fig. (4a) we demon¬ 
strate this through a sequence of time series of Sij for 
increasing noise strengths. We find that addition of weak 
noise delays the decay time of s^ gradually, eventually 
lowering its minimum to zero. This produces a delay in 
the onset of the next Sy peaks which in turn delays the 
input Ii to the post-synaptic neuron at a future time in¬ 
stant. The decay time of Sjj is increased further as the 
noise strength is increased. The periodically occurring 
maximal values of the inputs then arrive at the post- 
synaptic neurons later than in the previous cases and this 
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gets reflected in the firing pattern of neurons as depar¬ 
tures from the noiseless /previous values of the differences 
(ui — U 2 ) and those of the phase differences between the 
neurons. This can be seen more clearly in Fig.(4b) where 
the synaptic conductance for a single neuron is shown for 
increasing noise strengths. 

For a given set of p, rj, tr, the values of the in¬ 
stantaneous phases </>i and those of their temporal vari¬ 
ation, i.e., the instantaneous frequencies, are determined 
by the synaptic coupling strengths otijgij and the noise 
strength, which are inputs received at the synapse. When 
the variation of the instantaneous values with instanta¬ 
neous phases of the inputs to neuron 1 matches with that 
for neuron 2, then it is reasonable to expect CS to oc¬ 
cur in the EE system, because the instantaneous values 
(f> 1 and 4 >i of neuron 1 change in step with <f>2 and cj >2 
respectively of neuron 2. This forces the amplitudes of 
neurons 1 and 2 to become identically the same, because 
if they do not, then both conditions (f> 1 = $2 and </>i = fa 
cannot simultaneously be maintained. Such a plot (in¬ 
stantaneous frequency <f> versus instantaneous phase <j>) of 
the input to each neuron then has a characteristic flame 
shape (Fig.(5)), which exactly matches with that for the 
other neuron when CS occurs^, and this identity of the 
plots can be used as an indicator of CS. 


4.3 Input-Frequency curves and Lyapunov exponents 

The input-firing frequency curves for the coupled sys¬ 
tem plotted in Fig. (6) for different values of the coupling 
strength are typical characterizations of type-I neurons, 
differing from them only in that, even at Pi = 0.0 there is 
a non-zero frequency of neuronal firing, since Ii ^ 0. Lya¬ 
punov exponents have been calculated and shown along¬ 
side. 

These frequency-input curves indicate firing via a 
saddle-node bifurcation for the coupled EE system. 
Since the occurrence of a saddle-node bifurcation allows 
for a contraction region along and close to the stable 
manifold of the saddle, and in conformity with the result 
of2i, it seems reasonable to expect CS to occur among 
systems coupled through EE synapses. 

Fig. (6a) depicts identical neurons with EE coupling 
(dij = 9ji = 0.3 and other parameters as shown): the 
system is periodic with the presence of only a single fre¬ 
quency. Increasing the coupling strength of one neuron 
to gji = 0.45, all other parameters and conditions re¬ 
maining as before (corresponding to non-identical neu¬ 
rons with EE coupling) produces a second frequency af¬ 
ter Pi crosses a value of approximately 0.7 (Fig.(6b)). 
This is manifested in the Lyapunov exponent curve: the 
second largest Lyapunov exponent also becomes zero at 
the same point, indicating the presence of a second fre¬ 
quency in the system. As this transition from a periodic 
to a quasiperiodic state involves variation in and Pi, 
for given Ty , tr, a codimension two bifurcation is indi¬ 


cated at this point. 

4.4 Noise-induced variation of firing frequency 

Largest Lyapunov exponents in the presence of noise 
were calculated ini£ using the stochastic Runge-Kutta- 
4 method. We found that in both EE and IE cases, the 
largest Lyapunov exponent Ai becomes more negative on 
the addition of noise demonstrating the role of noise in 
bringing about order into the system. 

In the presence of weak noise, firing frequency changes 
with Pi differently for the excitatory and the inhibitory 
neurons, depending upon whether the coupling is of EE 
or IE type. For EE coupling, this is displayed in Fig. (7) 
and the plots for the noiseless case is shown in red in 
Fig. (7(b)). 

One observes that the frequency of firing of both ex¬ 
citatory neurons starts from a large non-zero value even 
at Pi = P 2 = 0 (Fig. (7)). As Pi is increased further, the 
frequency decreases rapidly, and after a certain value of 
Pi, the curve merges with that for <7 = 0, taking on the 
very same frequency values as for the noiseless case, so 
that from this point onwards, the frequency once again 
increases with increase of pi, but the increase is now grad¬ 
ual. 

4.5 Noise-induced synchronization & locking to partially 
synchronized state 

Numerous studies on CS in coupled systems have 
been previously reported in other systems (see for 
instance^— and references therein). Noise-induced 
CS in coupled theta neurons happens for EE coupling 
in a certain parameter range. However, we have not yet 
succeeded in observing it for IE coupling in a system of 
just two neurons. There appears to be definite relations 
between the magnitude of the output-difference between 
the two neurons averaged over all iterations (|ui — Zi 2 1), 
and the various parameters of the theory such as gij , Pi , 
and the noise strength tr, for both EE & IE couplings, 
and these may be quantified by plotting these quantities 
(Figs.(8)-(13)). 

In Fig.(8), (|ui — 7 / 2 I) is plotted as a function of gij for 
Pi = P2 = 0.1 in the case of EE coupling. It is seen 
that increasing the noise strength helps in bringing about 
synchronization for low values of gij . At higher coupling 
strengths, CS is destroyed as (|mi—W 2 I) increases steadily, 
but then it plateaus off, approaching a limiting value of 
approximately 0.5 at about g y = 6 for all noise strengths. 
Partial synchronization therefore takes place as the sys¬ 
tem gets locked to this state. 

Ini£ we found for EE synapses that for a given input, the 
critical coupling strength g c before which the neurons 
are completely synchronized, and beyond which they are 
desynchronized depends upon the noise strength as: 

g c ~ acr 1/2 (10) 
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where for the choice of parameters we had considered, we 
had a « 1.1. In that work, we also found for curves such 
as those in Fig.( 8 ), the following functional dependence 
of the synchronization error on g for the entire regime of 
coupling strengths after the onset of desynchronization: 

(|mi - u 2 \) ~ a{a) - b(a) - A (11) 

where the constants a and b are noise-dependent. 

The variation of (|ui — W 2 I) with pi for EE synapses 
can be seen in Fig.(9). CS does not happen for a = 0, 
but increasing noise strength decreases the difference. At 
lower nonzero noise strengths, the system oscillates be¬ 
tween the completely synchronized and the unsynchro¬ 
nized states. At a moderate noise strength of a = 0.6, 
the system shows CS for nearly the entire range of Pi 
shown in the figure. 

In Fig.(10), we show the variation of (|tti — U 2 I) with 
noise strength at different gij values, for pi = 0.1. We 
find that in neurons coupled through EE synapses, the 
transition to CS occurs at larger noise strengths for 
stronger coupling strengths; for weaker couplings, smaller 
noise strengths suffice for CS to occur. The system makes 
large oscillations between two metastable unsynchronized 
states before settling into the stable synchronized state. 
We find from an analysis of the plot in Fig.(10) that the 
point where desynchronization gives way to synchroniza¬ 
tion, i.e., (|iti — W 2 I) —> 0 begins to be approached, can 
be given by the relation 

Vturn » 0.5 ff 3 / 4 + P (12) 

We find that the maximum value that (|mi — W 2 I) takes is 
at minimal noise strength a and is related to the coupling 
strength g (for p = 0 . 1 ) by 

(pi-W 2 |)max~ 0.607 5 01 (13) 

This relation is made clear from the plot in fig.(10c). 

We now discuss our results for neurons coupled through 
IE synapses. 

5. IE SYNAPTIC COUPLING 

In Fig.(11a) we depict, as for the EE case, variations of 
both neuronal variables u\ and U 2 simultaneously, with 
P, g^ and r,j, in the absence of noise where the out of 
phase solutions are not the stable ones. The “firing” (ac¬ 
tive) points are more randomly distributed in this case. 
In contrast with the EE case, for IE coupling, we observe 
in Fig.(11a) that the distribution of the solutions about 
the diagonal u\ — 112 = 0 is far from being symmetric. 
Fig.(lib) shows the variation of Ui activity with P in 
a more standard way. The solutions for P less than the 
threshold value (ie for p < 0 ) show stable, non-oscillatory 
behaviour. All tt, for p > 0 however show large oscilla¬ 
tions. 


In the presence of noise, the IE system never achieves 
CS of its solutions. Introduction of weak noise does how¬ 
ever bring about some order even in the IE case and 
this can be visualised in the plots in Figs. (12) for noise 
strengths 0.5 and 1.0, where consecutive points in the 
U 1 -U 2 plane have been connected through a curve. The 
originally zig-zag trajectories in the noiseless case are now 
replaced by smooth curves which evolve in time explor¬ 
ing large regions of the phase space for increasing values 
of Pi in a definite, orderly manner. 

5.1 Analogy to a mechanical system 

In the case of IE synaptic couplings, addition of even 
very weak noise greatly affects the neuronal outputs, fir¬ 
ing rates and the inputs they receive. For this type of 
coupling, it is clear from the time series that the inputs 
have fast and slow branches always, characteristic of re¬ 
laxation dynamics. In this case, eqns.(5) and (6) can 
be rewritten (similar to the EE system) for N = 2, and 
a\ = 1 , 0:2 = — 1, as: 

«1 = /d(Ui) + Plfn{ui) -521(521 + l){u 2 )f n (ui) + £/„(«l) 
U 2 = /d(« 2 ) + p2fn{U2) + 012 (Sl2 + l){ui)fn(u 2 ) + £/n(«2) 
521 + h(u 2 )s 2 l + (<&(u 2 ) + ^(M2)5l2(5l2 + 1))S 2 1 
= -^(u 2 )S2lC(9 

Sl2 + h(ui)si2 + ($(«l) - Hi(lll)521 (521 + 1))S12 

= ->F(lli)si2£(9 (14) 

Here the nonlinear restoring force has variable spring co¬ 
efficient (<&(ui) — T(tti) 52 i(s 2 i + 1)). This restoring force 
is incremented for S 21 but reduces for S 12 . The system 
can be visualised crudely as two damped spring-mass sys¬ 
tems, one with a very stiff spring and the other with a 
very loose spring coupled together. The extension of one 
and simultaneous compression of the other spring leads 
to separation of time scales and relaxation oscillations of 
the synaptic conductances (Fig.(13)). These are fed as 
inputs to the neurons, again in the u 2 equation with a 
plus sign, while in the equation for u± with a minus sign. 
This discordance in the inputs results in non-synchrony 
of the neuronal outputs, and the excitatory neuron fires 
at a frequency different from the inhibitory neuron; this 
mutual difference in firing frequency is maintained even 
on increasing the common external constant input p: this 
is discussed further on in Section 5.3. Indeed, a canard 
solution (Fig.(13b)) appears to govern the underlying dy¬ 
namics. 


5.2 Variation of instantaneous frequency with 
instantaneous phase of neuronal inputs 

For IE coupling, though modulation of the firing 
thresholds by common noise occurs simultaneously for 
the two neurons, the mutual discordance of the inputs 
near the firing thresholds because of the opposite signs 
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in the respective coupling terms in eqns.(14) make it hard 
for the two neurons to cross the threshold and fire at the 
same points in time, as we discussed above, using the me¬ 
chanical analogy with the coupled damped spring-mass 
systems. Therefore although noise induces order even in 
IE systems (see Fig.(12)), yet complete synchronization 
becomes very hard to achieve because of the relaxation 
dynamics of the inputs: the “flame plots” for the two 
neurons do not match in this case (see Fig.(14) for ex¬ 
ample) . 

5.3 Input-Frequency curves and Lyapunov exponents 

Unlike the EE case, we cannot expect CS in the IE 
case, because the firing pattern (Fig.(15)) does not indi¬ 
cate presence of a saddle for both neurons — the firing 
frequency of the excitatory neuron shows non-monotonic 
behaviour with respect to the constant input /3: exhibit¬ 
ing very high firing frequencies at /3 ~ 0, then a sudden 
decrease at /3 ~ 0.3 and then a regular monotonic in¬ 
crease thereafter, and tending to maintain a nearly con¬ 
stant frequency difference with the inhibitory neuron for 
larger /3 in the deterministic system where a = 0. The 
nearly constant frequency difference however is sugges¬ 
tive of order emerging in their phase differences with in¬ 
creasing y. 

The plots in Fig.(15) are for parameter values identical 
with those in Fig. (6a) except for the coupling which is 
of IE type, (so that we now have now non-identical neu¬ 
rons): here again two different frequencies are present. 
In this case, for small /3j values, the large fluctuations in 
the frequency of the excitatory neuron (Fig.(15a)) when 
a = 0 arise because of the sign changes of the input 
{hi + YijLi a j9ji s ji(@j)) since ay = —1 ; the bifurca¬ 
tion parameter oscillates between the stable region where 
Ii < 0 and the unstable region /, > 0. Similar fluctua¬ 
tions are observed in the Lyapunov exponent curve also 
(Fig.(15b)), as shown in the inset curve, and the largest 
exponent even becomes positive in a small window of 
neuronal inputs. 


5.4 Noise-induced variation of firing frequency & delay of 
bifurcation 

For IE coupling, addition of noise causes the firing 
frequency of the inhibitory neuron to start from a 
high non-zero value just as in the EE case as the 
sign of the input term remains unchanged. However, 
for the excitatory neuron of Fig.(15), there is a non¬ 
monotonic dependence of its firing frequency on /3j (see 
Figs.(15),(16b)). In the absence of noise, when /3, = 0 
and aj = —1 then p < 0 and the neuron is below the 
firing threshold; hence the excitatory neuron does not 
fire. Introduction of weak noise however injects the 
required energy to the neuron to overcome the threshold 
and fire. In Fig.(16) we show that the frequency of 


firing in the presence of noise is higher than in the 
noiseless case (shown in red) but as the external input /3 
is increased, the effect of noise on the firing frequency 
dies down. 

We observe a very interesting phenomenon in the IE 
case when noise is introduced: the irregular firing of the 
excitatory neuron for small hi values which was accom¬ 
panied by two sharp peaks in the frequency (at the small 
hi values of 0.01 & 0.1) in the absence of noise, is re¬ 
placed in the presence of noise by a single broad peak at 
the much larger values of hi = 0-65 for er = 0.025, and 
hi = 1-6 for a = 0.075 (Fig.(16b)). Here noise appears to 
delay the bifurcation associated with this different firing 
behaviour. 

It is well known that the sweep of a control parameter 
which is slowly varying in time, causes a delay in the 
steady bifurcation points. The effect of Gaussian white 
noise into such a system reduces this delay in the bifurca¬ 
tion point depending upon whether the noise is additive 
or multiplicative^-. In eqns.(5) and (6) Gaussian white 
noise enters the system multiplicatively. The neuronal 
inputs especially for the IE synapse show relaxation os¬ 
cillations, having slow and fast time scales. This is quite 
distinct from the cases studied in~, where inputs vary 
very gradually and slowly with time. We note that the 
delay in the bifurcation associated with the change in 
the firing behaviour is for the excitatory neuron of the IE 
system, and not for the inhibitory neuron. 

As mentioned earlier, we observe from the time 
series that noise causes a delay in the synaptic decay 
tinted. While in the EE case this results in a change 
in the onset of the arrival of the maximal input values 
to the neurons and eventual synchronization (CS) 
of the neuronal outputs^ , in the IE case, because 
the signs of the input terms of the two neurons are 
mutually opposite, the neurons eventually are never 
able to fire in-step and hence CS is almost never achieved. 

Noise-induced delay of bifurcation in coupled theta 
neurons has been reported earlier in the literature^, but 
the coupling in their case was purely of EE type. 


5.5 Noise-induced synchronization & locking to partially 
synchronized state 

For IE coupling (Fig.(17)), the synchronization error 
(|ui — M 2 I) has a maximum magnitude at gtj =2.4 ap¬ 
proximately, for all values of noise strengths, and ap¬ 
proaches a saturating value of 0.27 at about <jij = 6 ir¬ 
respective of the noise strength, both for hi = 0 and for 
fii = 0.1. There are three interesting things to be noted: 
first, values of (|rti — U 2 I) are lower in the IE case in com¬ 
parison to EE coupling, but in the IE case too there is 
partial synchronization for large coupling strengths. The 
second point is that increasing noise strength beyond a 
critical value g c = 0.4 increases the difference between 
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neuron outputs while below this value the difference de¬ 
creases for nonzero cr, in contrast to the EE case where 
increasing noise strength brought down their difference. 
Finally, we see that for IE coupling in the noiseless sit¬ 
uation, the variation of (|ui — U 2 I) with 3 ^ at pi = 0.1 
is discontinuous at a small coupling strength of about 
gij = 0.3 (Fig.(17)). The reason for the discontinuity is 
not yet clear to us. 

Variation of (|ui — U2I) with the input /? for the IE case 
with n = r 2 , pi = P2 = P , 512 = 321 is shown in Fig.( 18 ) 
for different noise strengths. For the chosen values of 
parameters, we observe that there is a transition point 
at P = P c ~ 0 . 1 . For P > p c , (|ui — U2I) decreases 
with increasing a. There is a cross-over region between 
P = P Cl = 0.03 and p c , and for p < p Cl , (|ui — U2I) in¬ 
creases with increasing a. The plots show a dependence 
of the form: 

(|ui - u 2 \) ~ aP b , for p < p ci (15) 

where a & & vary with the noise strength a and 

(|iti - u 2 1) ~ kp + IP 1 / 2 + (2 - cr)/? 2 / 5 , for P > p c (16) 

where k , l vary from curve to curve and depend on the 
noise strength a. 

However, as mentioned earlier also, (|ui — M 2 I) never be¬ 
comes zero for IE synapses; instead, at large Pi values, 
the output difference saturates for all noise strengths to a 
limiting value of about 0.325 and the system gets locked 
on to this state. The discontinuity at pi = 0.1 for the 
noiseless case (in red) is a reflection of very small, non¬ 
negative values of the Lyapunov exponents (Fig.(15b), in¬ 
set) (discussed also in the previous section) because of the 
total input to the excitatory neuron oscillating between 
the stable regime when p < 0 and the unstable regime 
when Jj > 0 , depending upon the relative strengths of Pi 
and otjgjiSji, since ay = — 1 for neuron j. 

The variation of (|tti — U 2 I) with noise strength for IE 
synapses at different gij values, for /?j = 0.1 is displayed 
Fig.(19). Above a noise strength of about a = 0.15, 
increasing gij increases the difference between the neu¬ 
ron outputs; below cr = 0.15 also the system remains far 
from CS, but increasing g-ij can bring down the differ¬ 
ence. Comparing these IE plots with the ones in Fig.(10) 
for EE, we see that as usual EE and IE synapses show 
different behaviour. The uncoupled case gij = 0 is shown 
in both plots of Figs.(19),(10) for comparison. 


6. VARIATION OF SYNCHRONIZATION 
TIME WITH NOISE STRENGTH (EE 
COUPLING) 

In Fig. (20) we plot the time taken by the EE system 
to exhibit CS as a function of the noise strength, keeping 
all other parameters fixed. The normal clock function 
available with C-programming language has been used 


to measure the time T syn taken by the system to syn¬ 
chronize. 

To define T syn , the procedure we use is as follows. For 
a given value of the noise strength cr, the absolute value 
of the output difference \u\ — U2I is calculated, at each 
instance, for a large number k of iterations (k = 300000 ). 
Once |tti — U2I goes to zero in, say the itli iteration, 
(or more accurately, when this difference can be con¬ 
sidered to be negligible computationally - less than 
10 -6 ), the system is considered to have synchronized 
only if for all the consecutive and remaining (fc — i ) it¬ 
erations, |ui — U 2 I « 0 continues to hold good. When 
this synchronization condition is found to be satisfied, 
then the time (from the beginning of computation) at 
which |tii — U2I « 0 first occurred (at the ith iteration), 
is considered to be the synchronization time T syn . 

We observe that T syn is smaller for higher noise 
strengths. Infinitesimally small noise strengths suffice to 
completely synchronize two uncoupled neurons (g-ij = 0), 
though the time taken for this to happen may be longer. 
In this case, increasing noise strength causes CS to oc¬ 
cur for very small values of T syn . For increasing coupling 
strengths, there exist plateaus or bands of cr for which 
T syn is the same. For stronger couplings (c/12 = 321 = 0.5 
and 312 = 321 = 1-0) we observe the existence of 
metastable states in certain regimes of noise strengths 
where the times taken for CS to occur can be different. 
This suggests the possibility of coupled neurons utilizing 
the different synchronization times for encoding different 
features from the given inputs. Indeed, this opens up 
the possibility of exploiting “T syn switches” in artificial 
neural networks for feature extraction and recognition. 

7. CONCLUSION 

We have reported in this work results of some com¬ 
puter studies on noise-induced complete synchronization 
in systems of bidirectionally coupled type-I neurons. We 
have studied both EE and IE couplings. From analyses of 
these results, we obtain functional relations between the 
input strength, noise strength, coupling strengths and the 
distance from CS. Transitions from a completely synchro¬ 
nized regime to a partially synchronized state to which 
the system gets locked with increasing coupling strength, 
or with increasing input, other parameters kept constant, 
are observed. Similarly we observe that with increasing 
noise strength desynchronization yields to CS. We find 
that the critical noise strength where this occurs has a 
clear functional relationship to the coupling strength and 
the input. For IE coupling, we find that depending on 
whether the input /? is less or more than a critical value 
Pc the system’s departure from completely synchronized 
to a partially synchronized state occurs very differently. 

We find that Gaussian white noise can cause a delay 
in the bifurcation associated with the firing of the excita¬ 
tory neuron when the synapse is of IE type. We suggest 
that relaxation dynamics of the neuronal inputs in IE sys- 
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terns make CS of the outputs hard to achieve even in the 
presence of noise. Finally, we point out the possibility of 
utilizing the different times for complete synchronization 
as switches for extracting different features in networks 
of neurons. 
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Figure Captions 

Figure 1. (color online) Bidirectional coupling 
between two theta neurons. Left: E-E coupling, Right: 
I-E coupling. 

Figure 2. (color online) (a) Activity diagrams for EE 
coupling: Effect of different parameters on firing: (left: 
changing input fti = @2 = ft, middle: changing coupling 
strength g 12 = 321 = 9 , right: changing synaptic decay 
time T 12 = T 21 = t. 

(b) Plot in Fig. (2a, left) depicted in a more standard 
way, as variation of neuron outputs (on the y-axes) with 
the input ft (on the x-axis) showing bifurcation in the 
neuron activity. 

Figure 3. (color online) Complete synchronization in 
EE system in the presence of noise: left: cr = 0.375, 
right: cr = 1 . 0 . 

Figure 4. (color online) Noise produces delay in 
the decay time of the synaptic gating variables which 
eventually bring about CS in neurons with EE coupling. 
These plots illustrate this, (a) plots at the left are 
time-series of the synaptic conductances S 12 (in red) 
& s 2 i (blue) and the corresponding differences of the 
neuronal outputs (ui — u 2 ) are shown at the right for 
EE coupling. These plots show the approach to CS as 
cr is increased from top to bottom: a = 0 , cr = 0 . 20 , 
cr = 0.24, o- = 0.26, cr = 0.29002, a = 0.2900200338. 

(b) Detail of the plots on the left in (a) above are shown 
for increasing noise strengths (a = 0.0,0.20,0.24,0.26) 
for a single synaptic gating variable si 2 . 

Figure 5(color online) Instantaneous frequency vs 
instantaneous phase of inputs to neurons (“flame plots”) 
for EE coupling for: (a) cr = 0.0; (b) cr = 0.75. 
Parameters used are: fti = 0.0, /3 2 = 0.00001, 912 = 
521 = 0.3, t 12 = r21 = 2.0, tr = 0.1, 77 = 5.0, (Initial 
conditions: Q\ = 0.0, 9 2 = 0.0, si 2 = 0.0, s 2 i = 0.0). - 
Noise-induced CS occurs. 

Figure 6 (color online) Firing frequency and Lya¬ 
punov exponents versus input curves for EE coupling, 
for t 12 = t 2 1 = 2 . 0,77 = 5.0, t r = 0.1, cr = 0 : (a) 
.912 = 921 = 0.3 ; (b) g \ 2 = 0.3, 921 = 0.45. 

Figure 7 (color online) Frequency of firing & effect 
of noise on it for EE coupling for: (a) cr = 0.025 & 
cr = 0.05, (b) curves of the left graph plotted alongside 
values for cr = 0 . 

Figure 8 (color online) Variation of (|mi— u 2 |) with the 
coupling strength g \ 2 = g 2 \ = g for EE coupling. Param¬ 
eter values are: t \ 2 = r 2 i = 2 . 0 , fti = ft 2 = 0 . 1 , tr = 0 . 1 . 
(Initial conditions: 61 = 0.0, 0 2 = 0.01, S 12 = 0.0, s 2 i = 
0 . 0 ). 


Figure 9 (color online) Variation of (|tti — u 2 \) with 
the input ft for EE coupling: Parameter values are: 
T 12 = t 2 i = 2 . 0,912 = 921 = 0.3, t r = 0.1. (Initial 
conditions: 0 i = 0 . 0 , 0 2 = 0 . 01 , S 12 = 0 . 0 , S 21 = 0 . 0 ). 

Figure 10 (color online) (a) Variation of (|ui — u 2 \) 
with the noise-strength a for EE coupling, (b) crit¬ 
ical noise strength as a function of 9 (eqn.( 12 )); (c) 
(|iti - u 2 \) max versus 9 (eqn.(13)). 

Parameter values are: fti = ft 2 = 0.1, ti 2 = 
721 = 2.0, r r = 0.1. (Initial conditions: 

0i — 0.0, 0 2 = 0.01, si 2 — 0.0, s 2 i = 0.0). 

Figure 11. (color online) (a) Activity diagrams for IE 
coupling: Effect of different parameters on firing: (left: 
changing input fti = ft 2 — ft, middle: changing coupling 
strength g \ 2 = g 2 1 = g, right: changing synaptic decay 
time ti 2 = r 2 i = r. 

(b) Plot in Fig.(11a, left) depicted in the standard way, 
as variation of neuron outputs (on the y-axes) with 
the input ft (on the x-axis) showing bifurcation in the 
neuron activity. 

Figure 12. (color online) Noise-induced or¬ 
der in IE systems: for (from left to right): 

cr = 0.0, a = 0.5, a = 1.0. 

Figure 13. (color online) IE system: (a) relax¬ 
ation oscillations of inputs; (b) canard in the I\ - I 2 
plane; (c) time series of neuronal outputs iq. Here, 
912 = 921 = 0.3, 712 = 721 = 2.0, Tr = 0.1,7? = 5.0, fti = 
ft 2 = 0.1,0! = +1,02 = -I- 

Figure 14(color online) Instantaneous frequency vs 
instantaneous phase of inputs to neurons (“flame plots”) 
for IE coupling for: 

(a) ct = 0 . 0 , fti = ft 2 = 0 . 1 , 912 = 921 = 0 . 3 , 7 i 2 = r 21 = 
2.0, tr, = 0 . 1 , 7 ? = 5.0 (Initial conditions: Q\ = 0 . 0,02 = 
0.1,512 = 521 = 0.0); (b) (7 = 0.15, fti = ft 2 = 0.0, 
912 = 921 = 0 . 3 , 7-12 = 7-21 = 2.0, tr = 0 . 1 , 7 ? = 5.0, (Ini¬ 
tial conditions: 9i = 0.0, 9 2 = 0.01, S 12 = 0.0, S 21 = 0.0). 
CS does not occur for this case. 

Figure 15(color online) (a) Firing frequency and 

(b) Lyapunov exponents versus input curves for IE 
coupling for 712 = 121 = 2.0,?? = 5.0, tr = 0.1, a = 0, 
912 = 921 = 0.3. 

Figure 16 (color online) Frequency of firing as a 
function of constant input & effect of noise on it for IE 
coupling : for cr = 0, cr = 0.025, cr = 0.075. (a) inhibitory 
neuron, (b) excitatory neuron). 

Figure 17 (color online) Variation of (|iti — u 2 \) 
with the coupling strength gi 2 = g 2 i — g for 

IE coupling. Parameter values are: T 12 = 721 = 

2.0, fti = ft 2 — 0.1, tr = 0.1. (Initial conditions: 

0 i = 0.0, 02 — 0.01, S12 — 0.0, S21 — 0.0). 


11 




N. Malik, B. Ashok & J. Balakrishnan, Eur. Phys. J. B 74, 177-193 (2010). 
Noise-induced synchronization in bidirectionally coupled type-I neurons 


Figure 18 (color online) Variation of (|wi — W 2 I) with 
the input /3i = /?2 = P for IE coupling, (a) The curves 
correspond to eqn.(16) (/3 > /3 C ): k,l values associated 
with a given a are: Black: a = 0.0, k = 0.0564, l = 
-1.6731, Red: a = 0.3, k = 0.0281 ,1 = -1.3538, 
Green: a = 0.6, k = 0.0092,/ = —1.0569, Blue: 
a = 1.0, A; = —0.0052,/ = —0.6903 ; (b) Same as 
(a), but on a log-log plot, clearly showing cross-over 
behaviour. (c)Detail of plot on the left for /3 < (3 C ,curves 
correspond to eqn.(15) : a, 6 values associated with a 
given a are: Black: a = 0.0, a = 0.9820, b = 0.5042, Red: 
<7 = 0.3, a = 0.3305,6 = 0.1913, Green: a = 0.6, a = 
0.2150, 6 = 0.0481, Blue: cr = 1.0, a = 0.2215,6 = 0.0460. 
Parameter values are: T 12 = 121 = 2.0, <712 = 
(721 = 0.25, tr = 0.1. (Initial conditions: 
9\ — 0.0 ,62 = 0.01, S12 — 0.0, S21 = 0.0). 



p 

Fig. 2b (color online) 


Figure 19 (color online) Variation of (|rti — U 2 I) with 
the noise-strength a for IE coupling. Parameter values 
are: /3i = /?2 = 0.1, ti 2 = T 21 = 2.0, = 0.1. (Initial 

conditions: 9\ = 0.0, 9 2 = 0.01, S 12 = 0.0, S 21 = 0.0). 

Figure 20 (color online) Rate of complete synchro¬ 
nization T sync with noise keeping all other parameters 
fixed, for the EE system. 



Firing ■ > 


Firing > ■ 



M2 u, —> 

[ ■ 

-10 12 3 4 



0 0.2 0.4 0.6 0.8 1 


-10 12 3 4 


Fig. 3 (color online) 


(a) 


(b) 


Fig. 1 (color online) 


Filing ■ • ■ Filing ■ ■ ■ Firing ■ ■ ■ 



No Firing * * 

t 

No Firinq * • 







1 1 1 

0.8 

1 u 

■ ; 

0.8 

. 1 ■ ' '' 

'' 1 , 

0.8 



\ 





A 0.6 

' i ,s, 

a 0.6 

.. 1 * 

■>., ,, 

, . '9 

A 0.6 

1 

1 ' 

□ 0.4 

/ .4 

”0,4 

i'': * 

3 0.4 

1 ■ 


, / { 


1 1 1 1 »i 



0.2 

/. 

0.2 

' , ' "i 

0.2 

j 

1 , a 

H, 

0.2 0.4 0.6 0.8 

u,-> 

912 = 

0.2 0.4 0.6 0.8 

21 U,~> 


0.2 0.4 0.6 0.8 

u,~> 






101234 

-10 12 3 4 

0 5 

10 


Fig. 2a (color online) 


12 


















































































































SjjSa Sa.Sn $!;,% s ©s 21 


N. Malik, B. Ashok & J. Balakrishnan, Eur. Phys. J. B 74, 177-193 (2010). 
Noise-induced synchronization in bidirectionally coupled type-I neurons 


Fig. 4a (color online) 









O 4QO 1QCC) 1300 avn 2100 XOO 
Bmt— v 


13 

























































































































































































































































N. Malik, B. Ashok & J. Balakrishnan, Eur. Phys. J. B 74, 177-193 (2010). 
Noise-induced synchronization in bidirectionally coupled type-I neurons 




Fig. 6 (color online) 
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Fig. 15 (color online) 
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